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, We present a modified Thomas- Fermi theory that describes the increase of the hyperfine coupling 

constants of endohedrally enclosed atoms. We use the March boundary conditions corresponding to 
a positively charged spherical shell surrounding the nuclear potential to represent the effect of the 



fullerene shell. We obtain quantitative agreement with experimental data for N@Ceo and N@C7o, 
and find that fullerene radius dominates over the fullerene charge in its effect on the hyperfine 
coupling constants. We also present predictions for the hyperfine coupling constants of the endohe- 
dral nitrogen fullerenes between C60 and C500, and discuss the implications for proposed quantum 
computing schemes. 



i> 1 I. INTRODUCTION 

The nitrogen endohedral fullerene molecule, N@C6o, has received a great deal of experimen.t&)k2£££&l&& and 
theoreticaii2iii*i£ study. Unlike the metallofullerenes, the endohedral nitrogen acts essentially as a "free" atom with 



a5 = 3/2 quartet ground electronic state, although spatially restricted by a harmonic-oscillator- like potential within 
the fullerene Additionally, the fullerene acts a Faraday cage, shielding the spin of the endohedral nitrogen atom 
from stray fields^ and the sharp ESR spectra indicate long spin relaxation times^ which have inspired its use as a 

o : 

Q ' potential spin qubit for quantum computationii^iiSiiSiiLiSiiS Experimental^ and theoretical^ studies have indicated 



that suitable pulse sequences can reduce single qubit errors in these systems to the order of 10~ 6 , which is within the 
1CP 4 error threshold of Steanei^i 

An interesting physical phenomenon observed for endohedral atom systems is an increase in the hyperfine coupling 
constant as compared to its gas phase value, in particular, an increase of 54.1% for N@Cgo and of 49.1% for N@Grom^ 
To better understand the origin of this effect, Kobayashi et al. carried out ab initio MP2/uc-Huginaga+(2df) cal- 
culations, obtaining enhancements of 77% and 65% for N@Ceo and N@C7o, respectively^ In general, quantitative 
calculation of the hyperfine constant is a difficult problem, even for the free atom^ However, a clear qualitative 
model for the effects of changes in fullerene size and charge that allows prediction of these coupling constants would 
be helpful, not only for understanding the underlying physics, but also to provide estimates of potential sources of 
qubit error in applications to quantum information processing. In this paper we introduce a modified Thomas-Fermi 
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theory for the hyperfme coupling constants in endohedral fullerenes that is capable of quantitatively reproducing 
the experimental results and provides the desired predictive capabilities. We show that the model is accurate by 
reproducing the experimentally observed increase in coupling constants for N@Cgo and N@CVo . We then analyze the 
independent roles of fullerene cage size and number of electrons on this enhancement and make predictions for the 
hyperfme coupling enhancement for nitrogen atoms endohedral to fullerenes between Cqq and C§oo- 



II. THEORY 

Following the spherical boundary condition model for highly symmetric molecules developed by March^ and later 
applied to endofullerenes by Clougherty^ we treat the fullerene as a sphere of radius r = R, contributing n 7r-electrons 
and with a charge of Ne, around the nuclear point charge of Ze at the origin. For example, for N@C6o, Z = 7 and 
11 = 60. Transforming into the reduced coordinates, x, given by r — bx where 

1/9^ 1/3 



h= A\Jz) aBohr (1) 

and aBohr is the Bohr radius, we transform the cage radius R — > X. We then proceed to solve the usual Thomas-Fermi 
differential equation^ 

d 2 X X 3/2 



dx 2 x 1 ' 2 

with the additional boundary conditionj2i2ii 



(2) 



X '(X-)-x'(X+) = A (3) 

to account for the spherical shell, in addition to a continuity condition, = and the usual boundary 

conditions of x(0) = 1 and x(°°) — 0- Clougherty and co-workers^SI have examined the effect of treating the 
icosahedral (i.e., non-spherical) nature of the fullerene by means of a multipole expansion, but we neglect this effect 
in the current work. 

Next, we use the expression derived by Fermi , 28 ' 29 



dr (4) 



ri 

which relates the number of electrons, pi, having angular momentum Z, to the potential, V(r) — Zex(r)/r. Other 
expressions for determining the angular momentum assignments from Thomas-Fermi models have been developed, but 
give similar results for low Z. 30,31 The limits of integration are chosen so that the inteerand (and hence the integral") is 



a positive real number. In our case, we take the upper limit of integration to be some tunable cutoff radius, r c , inside 
the fullcrcnc shell, to count only the endohedral atom electrons and avoid contributions of the fullerene electrons to 
the integral. We are primarily interested in the number of unpaired I = electrons (vide infra). Considering only 
the "valence" I = electrons, there is no ambiguity due to Hund's rule, as there is for the I > electrons, so we may 
therefore express the number of fractional, i.e., non- integer, number of I — electrons as 



and the number of unpaired I = electrons as 



p{L7 = frac (^rr) (5) 



frac if frac < 1 



unpaired _ I Pl=0 ' 11 P 1=0 —2 , ,. 

Pl=0 - ) 1 _ J} frac . { frac 1 ■ W 

I, 1 Pl=0 ' 11 Pl=0 ' 2 

The contact-term of the isotropic hyperfine coupling constant 
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a = {—)g e gNHBLiNHS z ) l {*\p(p)\V), (7) 

is proportional to the electron spin-density (difference between spin up and spin down densities) at the nucleus, p(0) 
((S z ) ). However it is well known that the Thomas- Fermi method gives an infinite density at the nucleus, and 
qualitatively poor results near the nucleus (precluding extrapolation), so a direct evaluation of p(0) will failiSi A 
similar problem occurs in the treatment of the hyperfine coupling constants by semiempirical quantum chemistry 
methods, in which the 2s and 3s Slater orbital basis functions erroneously have zero density at the nucleus^ 

To avoid this problem, we present a heuristic argument in the language of atomic orbital theory. The atom-centered 
basis functions on a free atom may each be decomposed into a radial function multiplied by a spherical harmonic 
angular function. For orbitals with I > 0, a node occurs at the nucleus, resulting from the properties of the spherical 
harmonics, so the electrons in these orbitals have no density at the nucleus. This is unchanged when the atom is placed 
in a spherically symmetric potential, allowing us restrict our attention to the / = electrons only. The magnitude 
of the spin density at the origin, p(0), is the product of the probability density of finding an I = electron at the 
origin, times the number of unpaired I = electrons, p™™P aired ; given by Eq.©. Following the argument made by 
Pople et al.j& we can consider the basis, and hence the / = propability density term, to be unchanged by the 
chemical surroundings of the atom, and consider the changes to occur only in the basis orbital populations. Similar 
arguments have been made in the calculation of host medium effects (including endohedral fullerene inclusion) on the 
L/K electron-capture /3-decay ratio of 7 Bei^^ This allows evaluation of the ratio of p(0) for the free and endohedral 
atom cases, as in Eq.(|SJ, below, in which the unknown probability density term now cancels, leaving the ratio in terms 
of v ™ l P aired which we can calculate from Ea.Joll. 



Strictly, the fullerene carbon-atom basis functions may have a small non-zero value at the nitrogen atom nucleus. 
However, if we make a zero-differential overlap (ZDO) approximation, which is well justified by the lack of a chemical 
bond between the carbon and nitrogen atoms, this contribution is zero. In our calculation we neglect the contributions 
from fullerene electrons by considering large fullerenes, and by choosing the cut-off radius r c used to determine the 
integral in Eq.Q, to be sufficiently smaller than the cage radius, R, but sufficiently large to enclose the region one 
would chemically attribute to the endohedral atom. We discuss the choice of r c , with particular attention to the 
difference in applicability for N@C6o and P@Cgo in the next section. 

Following the above discussion, the hyperfine compression ratio of Buchachenko, 22 originally given in terms of the 
endohedral and free atom hyperfine coupling constants a and ao, respectively, may then be expressed as 

(8) 



unpaired/ j 7 j ?\ unpaired / r \ 

Pi=o {endohedral) — p [= ^ (free) 



(9) 



unpaired / r \ 

Pi=o (/ ree ) 

in which a consistent value of r c is used for the evaluation of both the endohedral and free atom terms, in order to 
make the probability density for the I = electron at the nucleus transferable between the two cases, as discussed 
above. 



III. RESULTS 

To determine the appropriate value of the cut-off radius, r c , we plot calculated values of <f> versus r c for N@Cgo 
and N@CVo in Figure^ obtained using the experimental fullerene dimensions given by Buchachenko. 22 To match the 
experimental measurement for N@C6o of <f> = 0.541, we take r c = 3.01 bohr. This yields values of (f> = 0.541, 0.562 and 
0.345 for Ceo, for 6.607 bohr radius (short axis) C70, and for 7.38 bohr radius (long axis) C70, respectively. Since the 
C70 molecule is prolate, the average value of <f> computed by weighting the short axis value twice is 0.490, in excellent 
agreement with the experimentally measured value of <f> = 0.491i2£ For comparison, we note that the values of <fi 
evaluated with r c — 3.00 bohr are 0.531 and 0.486 for Ceo and C70 (average value), and also in reasonable agreement 
with experiment. Our value of r c « 3 bohr is consistent with the spatial extent of the spin density distribution 
determined in the ab initio calculations of Kobayashi et ah 12 However, our attempt to perform a similar calculation 
for P@C6o was unsuccessful, due to the larger size of the P atom. This requires a larger r c and thus unavoidably 
introduces possible contributions from fullerene electrons. Thus we restrict the analysis below to nitrogen endohedral 
fullerenes. 
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We have examined the relative roles of the number of cage 7r-electrons, n, and radius, R, on cj>, shown in Figure El 
Intuitively, one expects <f> to be inversely proportional to fullerene radius R (at constant number of cage electrons n), 
and to be directly proportional to n (at constant R). Variations in n and R may arise experimentally by endohedral 
inclusion into defect fullerenes^ by use of boron or nitrogen substitutions in place of carbon atoms in the fullerene, 37 ' 38 
or as the result of interactions with surfaces^2a£ In practice, simultaneous variations of n and R are unavoidable. 
Nevertheless, we will analyze ideal independent variations of these two parameters in order to establish their relative 
importance. The papers by Buchachenkc 22 and by Kobayashi et al.fH both make the qualitative statement that 
cage radius (R) is more important than the cage electrons (n) in enhancing the hyperfine coupling of the endohedral 
nitrogen atom. The simplicity of the present model allows us to separate and directly test these two contributions. 
Figure [21 confirms these statements, as is visible from the much steeper slope of the constant-n curves (i.e., varying 
R), as compared to the constant-i? (i.e., varying n) curves. For N@C6o, at An = and AR = 0, the slopes are 
d(f>/dR = —0.35 bohr -1 and d(f>/dn = 0.0031. We discuss the implications of the magnitude of these dependencies for 
quantum computation in the next section. 

To examine the hypothetical endohedral fullerenes N@C„ with 60 < n < 500, most of which have not yet been 
systematically studied, we have assumed that the ratio of n/S is constant, where S = 4irR 2 is the surface area of a 
spherical fullerene. Thus for C 60 , S = 543 bohr 2 and n/S = 0.11; for C 70 , S = 592 bohr 2 and n/S = 0.12. Figure 
[3]sliows the behavior of the hyperfine compression ratio, </>, for N@C„ with 60 < n < 500, for n/S — 0.10, 0.11, and 
0.12. As expected, (j> is seen to be a monotonically decreasing function of n, regardless of the choice of S/n, and 
asymptotically approaches zero (i.e., the free-atom limit of no enhancement of the hyperfine coupling constant) as the 
fullerene becomes larger. Depending on the exact value of S/n, we find (j> < 0.1 for 130 < n < 160 and <j> < 0.01 for 
350 < n < 380. Even for the hypothetical C500 fullerene, with a radius of ~ 40 bohr, we find a finite enhancement of 
0.002 < 4> < 0.004, which could be detectable experimentally. 22 

IV. SUMMARY AND DISCUSSION 

Using a modified Thomas- Fermi theory, we have developed a simple physical model for the effect of endohedral 
inclusion on the hyperfine coupling constant of nitrogen. The only free parameter in our model is the cutoff radius, 
r c , used to partition the atom from the fullerene. A single choice of r c « 3.0 bohr is found to give quantitative 
agreement with recent experimental results for both N@C6o and N@C7o- Moreover, the simplicity of our model allows 
us to separate the relative contributions of cage potential and size on the hyperfine coupling constants, and to thereby 
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identify the fullerene radius as the dominating factor in the scaling of the hyperfine compression ratio cj>. We have 
also predicted the hyperfine coupling constants for nitrogen endohedral to hypothetical fullerenes as large as C500 and 
find that a finite enhancement should exist up to this size. 

Spin exchange betwen the electron and nuclear degrees of freedom mediated by the hyperfine interaction is necessary 
to implement the quantum cellular automata scheme of Twamley^i^ Since global operations on all the qubits are used 
to evolve the computational unitary in this scheme, deviations in the hyperfine coupling constants of the individual 
sites, e.g. due to deformations resulting from fullerene interaction with the substrate surface, would require additional 
complication in the nuclear-electronic CNOT and swap operations in order to be robust against these variations. Our 
results in Section IIIII indicate that the hyperfine coupling is quite sensitive to these deformations, and that even a 
1% (0.07A) change in the N@Cgo radius would result in a ~ 5% change in the compression ratio 0, corresponding 
to an approximately 0.5 MHz shift in the hyperfine constant. Using the fidelity measure F = |Tr(V[/t)|/Tr(J7[/'f), 
where U is the desired unitary and V is the actual (erroneous) unitary^ a 0.5 MHz shift away from the nominal 
hyperfine interaction gives 1 — F = 2 x 10~ 4 for the electron-nuclear CNOT operation. Atomistic calculations treating 
the fullerene-surface interactions in more detail, may provide useful estimates for the extent of this deformation on 
experimentally relevant surfaces. 
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FIG. 1: Hyperfine compression ratio, <f>, as a function of cutoff radius, r c . For r c — 3.01 bohr, the compression ratios calcualated 
for N@C6o and the appropriate weighted average for N@C7o are 0.541 and 0.490, respectively, in excellent agreement with the 
experimental values 0.541 for N@C R o and 0.491 for N@C7tn^ 
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FIG. 2: Effects of changing n (positive slope curves) and R (negative slope curves) for Ceo (no = 60, Ro = 6.578 bohr) and a 
spherical model C70 (no = 70, Ro = 6.8647 bohr) on the hyperfine compression ratio (j>. Here n — no + An and R — Ro + AR, 
where no, -Ro are the equilibrium values of n and R, respectively. For Ceo at An = and AR = 0, dcj>/dn = 0.0031 and 
d<j>/dR = -0.35 bohr" 1 . 
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FIG. 3: Values of the hyperfine compression ratio, cj> for the N@C„ with 60 < n < 500, assuming surface charge/area ratios, 
n/S, of 0.10, 0.11, and 0.12. 



